function QuickXglomPipeline(img,MRorCT,ThresholdRange,PixelWidth,DatasetID)

[stored_hist stored_diams ThresholdRange diams_range] = ThresholdingProcess(img,1,ThresholdRange);
GlomStats = AnalyseThresholdStatistics(stored_diams, stored_hist, PixelWidth, ThresholdRange, diams_range,0)
%PlotXGlomOutputDist(stored_hist, stored_diams,DatasetID,PixelWidth,ThresholdRange);
%
%PlotXGlomOutputDist_version2(stored_hist,stored_diams,DatasetID,PixelWidth,ThresholdRange, diams_range,0,GlomStats)
PlotXGlomOutputDist_version3(stored_hist,stored_diams,DatasetID,PixelWidth,ThresholdRange, diams_range,0,GlomStats)
